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Abstract: We provide a theory for employing Forster resonance energy transfer (FRET) 
measurements to determine altered heteropentameric ion channel stoichiometries in 
intracellular compartments of living cells. We simulate FRET within nicotinic receptors 
(nAChRs) whose a4 and P2 subunits contain acceptor and donor fluorescent protein 
moieties, respectively, within the cytoplasmic loops. We predict FRET and normalized 
FRET (NFRET) for the two predominant stoichiometries, (a4)3(|32)2 vs. (a4)2(P2)3. 
Studying the ratio between FRET or NFRET for the two stoichiometries, minimizes 
distortions due to various photophysical uncertainties. Within a range of assumptions 
concerning the distance between fluorophores, deviations from plane pentameric geometry, 
and other asymmetries, the predicted FRET and NFRET for (a4)3(P2)2 exceeds that of 
(a4)2(P2)3. The simulations account for published data on transfected Neuro2a cells in 
which a4p2 stoichiometries were manipulated by varying fluorescent subunit cDNA ratios: 
NFRET decreased monotonically from (a4)3([32)2 stoichiometry to mostly (a4)2(P2)3. The 
simulations also account for previous macroscopic and single-channel observations that 
pharmacological chaperoning by nicotine and cytisine increase the (a4)2(P2)3 and 
(a4)3(P2)2 populations, respectively. We also analyze sources of variability. NFRET-based 
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monitoring of changes in subunit stoichiometry can contribute usefully to studies on 
Cys-loop receptors. 

Keywords: nicotine; cytisine; NFRET; nicotine addiction; Parkinson's disease; ion channels 

Abbreviations: EGFP: enhanced green fluorescent protein; ER: endoplasmic reticulum; 
FRET: Forster resonance energy transfer; nAChR: neuronal nicotinic acetylcholine 
receptors; mCherry: monomeric cherry fluorescent protein; NFRET: normalized Forster 
resonance energy transfer 



1. Introduction 

A superfamily of ligand-gated ion channels comprises homo- and heteropentameric combinations of 
subunits. The superfamily includes neuronal nicotinic acetylcholine receptors (nAChRs) comprised of 
a (a2 to alO) and P (|32 to [34) subunits. Possible changes in the stoichiometry of these subunits within 
the a4p2 pentamer have generated interest, because the (a4) 2 (P2) 3 and (a4)3(P2) 2 stoichiometries differ 
in their sensitivity to agonists [1,2] and to antagonists [2], their rectification properties [3], their Ca 
permeability [4], their sensitivity to upregulation by chronic nicotine and other drugs [5,6], and 
possibly their subcellular localization. Some of these characteristics may provide signatures for 
determining the stoichiometry of nAChRs when they appear on the plasma membrane. However, 
several nAChRs and other Cys-loop receptors localize to the endoplasmic reticulum (ER) [6-8], where 
stoichiometry measurements with biochemical methods are either tedious [5] or impossible. Therefore 
it is necessary to have additional robust measurements of changes in intracellular receptor stoichiometry. 

Chronic nicotine increases the assembly of (a4) 2 (P2)3 receptor pentamers within the ER [6,9], 
enhances the ER exit of nAChRs, and therefore upregulates plasma membrane receptors [6]. Such 
selective intracellular pharmacological chaperoning of acetylcholine receptor number and stoichiometry 
explains several aspects of nAChR upregulation by nicotine and other nicotinic ligands [6,9-15]. 
In addition to the obvious connection with nicotine addiction, pharmacological chaperoning may 
provide the mechanistic basis for the inverse correlation between a person's history of smoking and 
his/her susceptibility to Parkinson's disease [12]. These points increase our interest in intracellular 
measurements of alterations in a4[32 nAChR stoichiometry. 

We and others have utilized fluorescent protein (FP)-tagged a4 and P2 nAChR subunits and 
fluorescence microscopy in transfected cells to study altered a4p2 stoichiometry [6,16-19]. This paper 
provides the theoretical basis of an improved FRET method for monitoring these changes. In our 
previous experiments, we incorporated donors in 50% of the a subunits, and acceptors in the 50% of 
a subunits, leaving the P subunits unlabelled, or vice-versa. In our more recent experiments, 
all P subunits contained donors, and all a subunits contain acceptors, or vice-versa [6]. In this case, (1) 
in a pentamer with three acceptors and two donors, FRET is stronger than in the two acceptor-three 
donor case; and (2) FRET is also stronger between adjacent than between non-adjacent subunits. The 
key prediction is the ratio of the FRET for the stoichiometries, (a4)3(P2) 2 vs. (a4) 2 (P2)3. We show that 
the new method does not require additional strong assumptions about symmetry, or about distances 
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between fluorophores. The method may be considered a special case within the tradition of analyzing 
pentameric proteins by FRET [20-23]. The theory is of interest because the measurements have 
yielded estimates about changes in stoichiometry that are consistent with known structural information 
and with chaperoning by nicotine, as well as a novel effect of chaperoning by cytisine. 

2. Results and Discussion 

2.1. Simulations 

We calculate the relative FRET efficiency (E) for nAChRs composed of the (a4-mCherry)3([32-EGFP)2 
vs. (a4-mCherry)2([32-EGFP)3 stoichiometry. We define the average FRET values as 2 and E 2 3 

respectively, so that the desired ratio is E i2 l E 2i . We then extend these calculations to NFRET and 

calculate the analogous ratios, NFRET^^INFRETj^. The assumptions provide a special case of 
previous theories about pentameric proteins, but are less restrictive than previous analyses of 
pentameric ion channels [18,21], as described below (see Figure 1). 

(a) All [32-EGFP subunits have identical structure, and all a4-mCherry subunits have identical 
structure. For simplicity, we generally omit the description of the fluorophore when discussing 
stoichiometry; thus in most cases (a4) 2 (P2) 3 implies (a4-mCherry) 2 (P2-EGFP) 3 . 

(b) FRET interactions are fully determined by relative positions of subunits within the pentamer. 
The EGFP donor fluorophores need not be located at the same radial distance as the mCherry 
acceptors. The donor fluorophores also need not have the same dipole moment angle with the radius, 
nor lie in the same plane, as the acceptors. Indeed, the intracellular domain of the wild-type a4 subunit 
has roughly twice as many amino acids as that of the wild-type [32 subunit, rendering it likely that the 
a4-linked fluorophores have a different disposition from [32-linked fluorophores. 

(c) Donor subunits can be either adjacent or non-adjacent to acceptor subunits. Figure 1 A depicts a 
receptor with (a4)2(P2) 3 and (a4)3((32) 2 stoichiometry in the upper and lower panels respectively, and 
also shows the nomenclature for distinguishing among the three or two donor fluorophores, so that 
calculations can proceed for each donor. Figure 1A also shows the usual assumption, that two 
interfaces exist in which the [3 subunit is immediately counterclockwise to the a subunit. These 
appropriately polarized interfaces, which we term [32 < a4, are important for function because they 
form high- affinity agonist/antagonist binding sites; however the presence of a ligand site does not 
affect the FRET properties (the analogous terminology applied to GABA A or GluCl receptors would 
be a<[3, because for such receptors, the principal interface for binding lies on the [3 subunit). The 
fluorophores are shown as green (EGFP donor) and red (mCherry acceptor) bars, to emphasize the 
importance of the fluorophores' dipoles. The locations and orientations of the dipoles are arbitrary. 
The wild- type a4 M3-M4 intracellular loop has a sequence roughly twice as long the wild-type [32 
subunit; this loop is depicted schematically, and the EGFP donor fluorophore projects from the loop. 
Figure 1A thus emphasizes that the donor and acceptor fluorophores can also differ with respect to 
their radial distance. The donors can also reside in a different plane from the acceptors; this point is not 
shown. The distances between donor and acceptor fluorophores, angle between the donor emission and 
acceptor absorption dipoles, and spectral overlap integrals, lead to four possible microscopic FRET 
probabilities pj for each donor. These are calculated as described below. Figure IB depicts the energy 
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transfer pathways when a donor is flanked by adjacent acceptors; and Figure 1C shows the pathways 
when two acceptors are both non-adjacent to a donor. Depending on the order and stoichiometry of the 
subunits, only specific subsets of these pathways exist for each donor. Figure ID tabulates the 
"^-factors" for the nAChR stoichiometries under consideration. The subunit order is clockwise, seen 
from the extracellular surface. The two or three [32-EGFP subunits are numbered, as in A. The first 
column corresponds to A, top panel; the second, to A, bottom panel. Each cell in the table reports «y 
values that equal 1 (all others = 0), corresponding to possible energy transfer pathways. The «y values 
control the appearance of terms in Equations 2 through 6. 



Figure 1. Diagrams depicting FRET analysis of nAChR stoichiometry. (A) Two examples 
of stoichiometry, thought to correspond to the major a4|32 nAChRs in brain [12]. Both 
stoichiometries contain two [32 < a4 interfaces. The < character implies a correctly 
polarized high-affinity ligand binding interface, contrasting with the arrows in (B) and (C) 
which denote energy transfer. The nAChR is viewed from the extracellular surface. In the 
top panel the stoichiometry is (a4) 2 (P2) 3 , and in the bottom panel (a4) 3 ((32) 2 . The donor 
molecules are labeled 1 through 3, as in the equations in the text; (B) Nomenclature for 
energy transfer to acceptors that are adjacent to the donor. The first subscript denotes the 
identity of the [32-EGFP donor subunit. The second subscript is ±1 for adjacent and 
non-adjacent energy transfer, respectively; the positive sign denotes clockwise transfer. 
[32-EGFP donor #1 is shown with two immediately flanking a4-mCherry acceptors. 
Because these two pathways exist, the "^-factors" n\ t -\ = = 1. The arrows show the 
fluorophore separations R±\, R±2 and the corresponding energy transfer probabilities, p±\, 
p±2 Clockwise energy transfer is shown as a black arrow; anticlockwise, as gray. The 
remaining two subunits are shown in gray; they may be either donors or acceptors; 
(C) Nomenclature for energy transfer to acceptors that are non- adjacent to the donor. 
[32-EGFP donor #1 is shown with two a4-mCherry acceptors that are non-adjacent to the 
donor (but incidentally adjacent to each other). Because these two pathways exist, the 
"n-factors" »i,-2 = Wi,+2 = 1. Other details as in (B); (D) Table showing the "n-factors" for 
the nAChR stoichiometries under consideration. See text. 

A 

P<app<a 




1 23 12 
(3<app<a p<ap<aa 

n-factors" that equal 1 (all others = 0) 
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For Equations (1) through 8 below, we consider only those donor and acceptor molecules within 
assembled receptors. The rate constant t d \ , for decay from the excited state increases linearly with the 

contribution from each potential donor-acceptor interaction: 

_ _ R 6 

T'daj = *d (1 + fp E n i.jPj ) = t d (1 + /pE Hj (1) 

K j 

where is the unperturbed rate constant. The equations select among the available energy transfer 

pathways as follows: each pj is multiplied by a parameter that we term an "^-factor" mj. The index i 
describes one of the two (i = 1, 2) or three (i = 1, 2, 3) individual donor fluorophores in the pentamer, 
and j = -2, -1, +1, +2 describe acceptors located at adjacent (±1), non-adjacent (±2) positions, where 
negative numbers indicate anticlockwise transfer as described in Figure 1B,C. The "n-factor" ntj is 
unity if donor i has an acceptor at relative position j, and zero if not. For example, if there is an 
a4-mCherry acceptor immediately clockwise to [32-EGFP donor #2, then «2,+i = 1; and the associated 
p+\ therefore exists in the term describing [32-EGFP donor #2 in Equations 2, 3, 9, and 10 below. 
Figure ID is a table of all "^-factors" that equal unity. 

Measured spectral FRET values always return products of the Forster efficiency E multiplied by 
the fraction fp of fluorophores participating in the energy transfer within assembled pentamers. 
The parameter is less than unity because donor and/or acceptor fluorophores fold incompletely, 
mature incompletely, become proteolyzed, undergo partial bleaching, or otherwise fail to participate in 
FRET [24]. The parameter fp is included in Equations 2 and 3 below. The possibility of unpaired, free 
donors residing outside of pentamers is considered later. 

In Equation 1, the FRET probabilities are calculated, as usual, by terms of the form p . = — . Rja 

R j 

is the Forster distance for each y'-valued donor-acceptor pair, and Rj is the distance between each 
y-valued donor-acceptor pair (length of the arrow in Figure 1). The Forster distance is determined by 
the spectral overlap between EGFP and mCherry (this is constant among all y'-valued pairs) and by the 
dipole angle between donor and acceptor; the latter parameter determines the orientation factor k , 
which potentially varies among the j values. For the specific assumption that k = 2/3, the Forster 
radius Ro for the EGFP-Cherry pair is 51.4 A [25]. However, we do not explicitly assume that the 
fluorophores move so quickly that they sample many orientations; indeed, this assumption is 
unacceptable for fluorescent proteins embedded within a loop. The parameter A, described below, 
embodies the possibility that k differs among donor-acceptor pairs within a pentamer. 

We now compute the desired ratio, £3, 2/E2, 3- The (a4) 2 (P2) 3 nAChR stoichiometry of Figure 1 A, 
top panel, has 3 donor molecules: 

= T D [l + fp (P + l + P-l )\ *DA,2 = *d [(1 + f P (P + l + P-2 )1 

?DA,3 = ^ + fp(P-l+P + l)] ^ 

In the (a4) 3 (P2) 2 configuration of Figure 1A, lower panel, the two donors have 

^DA,l ~ t D %+fp{p + 1 + P-1 + P+2 )l> 1 ~DA,2 ~ 1 ' D [l + /p(p+i +P-1 +P- 2 )] (3) 

We now apply the equation, 
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,-, , 7 DA,i 

E,=l-— (4) 

1 D 



We average over the two donors or the three donors, as appropriate (all conventional optical methods 
perform such averaging, because the resolution is much poorer than the subunit spacing). Thus, 



£■(3,2) = !-— ^- = 1 ' 



T D 2 



1 1 

■ + ■ 



1 + f P (P + l + P-l + P + l) l + f P (P+l + P-l + P-2 ). 



(5) 



and 

£(2,3) = 1-| 



- + -, ^ + 



l + f P {p + i+P-i) l + fp(p + i+P-i) ^ + f P (P- l +P +2 ) 



(6) 



These relations allow us to calculate E 32 /E 23 . 

(d) An unknown value is the distance between fluorophores on adjacent subunits, because there are 
no structural data for the cytoplasmic M3-M4 loops. We define an average distance R\^ ve between 
fluorophores on adjacent subunits. We then calculate the quantities of interest for a range of R\ fl ve 
values. The sensitivity of FRET measurements is limited to distances within a factor of two of Rq. The 
lower limit of our simulations, 25 A, represents both ~ Roll and the closest possible distance between 
adjacent fluorophores of the GFP family. The upper limit corresponds to ~2R 0 = 100 A. As discussed 
below (Figure 2), the E i2 1 E 2J and NFRETijJNFRETj?, ratios both approach the theoretical limit of 

1 .5 at -80 A, obviating exploration of greater R\, ave values. 

In order to perform specific calculations in the absence of structural data about the M3-M4 loops, 
we further define a regular pentagon-shaped receptor, as well as variations around this standard, 
as follows. 

(e) If the five fluorophores are located at the vertices of a pentagon, we may write, 

^+1 = ^-1 = R l,ave > ^+2 = ^-2 = R 2,ave = GR \,ave (7) 

We define a "geometry factor", G, which equals 1 .62 (the ratio of the length of a diagonal to a side 
of a pentagon) for the pentagonal receptor. Any departure from a regular plane pentagonal structure, or 
any difference between the plane of the a4-mCherry vs. [32-EGFP fluorophores [26], would decrease 
G. Figure 2B shows how the calculated E i2 l E 23 values vary for G between 1 .3 and 1 .62. 

(f) We define an "asymmetry factor", A, that can arise from unequal Rj values, and/or of k values, 
between the donor and two possible flanking acceptors that are both adjacent to the donor 
(as in Figure IB). For simplicity, we allow the same parameter to represent the asymmetry for the two 
possible acceptors that are non-adjacent to the donor (as in Figure 1C). 

P +l = (1 + A)A, flve . P-l = (1 - A )Pl,ave , P + 2 = (1 + A )Pl, aV e . P-2 = (1 - A )P 2 ,ave 
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Figure 2. Simulations. (A) The j-axis is the simulated E i2 / E 2J , the ratio of classical 

FRET efficiency for the two stoichiometries. The vertical arrows give parameters 
associated with the regular pentagonal receptor. The theoretical curves are computed for 
various values of the average distance between fluorophores on adjacent subunits Ri, a ve- 
X-axes are (left) the geometry factor G, and (right) the asymmetry factor A. The assumed 
fp = 0.5; (B) Same parameters as in (A); the plotted values are NFRET^INFRET2^, 
the ratio of NFRET efficiency for the two stoichiometries; (C) (D) Values of 
NFRET X2 /NFRET 2 ,3 vs. asymmetry factor A for assumed//. = 0.2 (C) or 0.8 (D). 
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The special case of A = 0 defines a regular pentameric, pentagonal receptor: (i) the five 
fluorophores are located at 72° intervals with respect to the nAChR axis of pseudo symmetry, 
presumably the channel pore; and (ii) their dipole moments make a common angle with the radius. 
In the simulations, we allow the range A = 0 to 0.6 (the latter value indicates that p +l differs from 
p+\lp~\ ~3 p-i by 4-fold). Figure IB shows an extreme case in which the difference between R-\ and 
R+i produces p+\lp-\ slightly greater than the 4-fold range. Monte Carlo calculations on a regular 
pentagon (Figure 10C of Reference [21]) show that k between two adjacent subunits lies mostly 
within a range of ±10%; the k values between non-adjacent subunits are somewhat larger, but again 
lie mostly within a range of ±10%. The corresponding range is p+\lp-\ ~3. 

These relations lead to Figure 2 A, which plots E 32 /E 23 versus G and A; this curve is repeated at 

5 A intervals of R\, a ve- Our experiments have used sensitized emission measurements, in which FRET 
is assessed by acceptor fluorescence [6,18,27]. The average total fluorescence from a given pentamer 
is the average of the energy transferred via resonance, times the quantum yield of the acceptor 
fluorophores Qa(. 

1 daj Ma 1 

*D 

Thus, the Qa% factor cancels out, and the equations for sensitized emission FRET are the right-hand 
sides of Equations 5 and 6. We also define f D , the fraction of donors that actually reside in fully 

assembled pentamers (including fully assembled, but immature, non-glycosylated or post-translationally 
modified pentamers). The remaining donors would reside in FRET-incapable soluble proteins or 
in partially FRET-capable, partially assembled receptors. The prime notation reminds us that 
FRET capability is additionally represented, within assembled pentamers, by fp. Assuming that f D 

does not depend on stoichiometry, this factor does not contribute to the FRET ratio calculations for 
sensitized emission. 



2.1.1. Calculations of NFRET 



Most of our experiments use NFRET, a variant in which the sensitized emission is normalized 
as follows: 

NFRET = J^ A , /o\ 




In Equation 8, Ia is the fluorescence from all acceptor molecules, both those that receive energy 
from donors and those that do not [24]. To the extent that I A is measured accurately by direct excitation 
at the excitation wavelength for the acceptor, I A does not depend explicitly on the stoichiometry and 
therefore cancels out when one takes the ratio NFRET^^INFRETj^. In Supplementary Material, we 
analyze possible deviations from this assumption. 

In Equation 8, Id is the donor fluorescence. This fluorescence is reduced because some EGFP 
molecules undergo FRET within pentamers; but EGFP molecules outside pentamers have much less or 
zero FRET. Therefore we write, 
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h=c^-f D E) (9) 

where c is a multiplicative constant which cancels out when ratios are computed. Id decreases to zero 
if both E — 1 and f D = 1 . Our simulations use f D = 1 , because this is the most "pessimistic" 

assumption from the viewpoint that NFRET is a reasonable monitor of FRET. With these assumptions, 
NFRET is approximately proportional to FRET efficiency E for E < -0.75, but NFRET becomes 
infinite as E approaches unity [28]. This poorly behaved range is not approached in our simulations, 
for three reasons. First, f P remains <1. Thus, in our previous studies, the maximal FRET efficiency was 
48% [17,18,29]. Second, R lave remains > R ( /2. Third, the process of computing NFRET X2 INFRET 2 ^ 
cancels out additional poorly behaved characteristics of NFRET. As a result, Figure 2 reveals that 
NFRET32/NFRET2J varies in much the same way as E i2 1 E 2i over a wide range of simulated 

parameters, as long as fp < -0.8 (NFRET^^NFRET^ does become poorly behaved for 
fp > ~0.9). Our techniques do not require that NFRET remains strictly proportional to FRET over the 
simulations, or that NFRET 32 1 NFRET 2i =£ 3 2 /E 23 , but simply that NFRET 3a /NFRET 2 j lies in a 

useful range for all plausible values of nAChR structure. 

Measurements of net FRET are treated in Supplementary Data. We found measurements of net 
FRET less useful those of NFRET. 

2.1.2. Summary of the Predictions 

We emphasize the major predictions that 1 < E 3 2 I E 2 3 < 1 .5 and that 1 < NFRET^^/NFRET^ < 1 .5 for 

all reasonable structures of a4[32 nAChRs. The limit of 1.0 at small values of Ri^ve occurs because the 
energy transfer approaches equal efficiency for all donors, regardless of the number of acceptors. The 
limit of 1.5 at large R\, arises from the fact that the average [32-EGFP donor has 1.5 times as many 
adjacent a4-mCherry acceptors in the (a4)3([32)2 stoichiometry as in the (a4) 2 (^2)3 stoichiometry; and 
at distances much greater than R 0 , the adjacent subunits dominate the energy transfer (p±\ □ p±2). The 
quantitative predictions are rather more sensitive to variations in the asymmetry factor A than in the 
geometry factor G (Figure 2). 

2.2. Agreement with Biased Transfection Experiments 

We have utilized nAChR a4-mCherry (acceptor) and [32-EGFP (donor) subunits. These 
fluorophores have, in our hands, optimal characteristics for pixel-by-pixel sensitized emission in living 
cells. These proteins are expressed in Neuro-2a cells, a favorable system in which membrane protein 
expression remains linear with respect to several factors under the experimenter's control [6,16-19]. 
In one set of experiments, Neuro-2a cells were transfected with 0.33, 0.5, or 0.67 mole fraction cDNAs 
for the [32-EGFP subunit (the balance was a4-mCherry cDNA, for a total of 1 ug plasmid DNA) [27]. 
These give more robust signals than in our previous studies utilizing only fluorescent a subunits, or 
only fluorescent [3 subunits. Measurements based on individual pixels (rather than on the image from 
an entire cell) make the best use of available data from confocal imaging [30-32]. For instance, despite 
the substantial pixel-to-pixel variance, the mean NFRET values have standards errors of the mean 
much less than 1%. Sensitized emission measurements are nondestructive and can therefore be repeated 
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many times on a single cell. NFRET equals Ida normalized to the square root of the donor x acceptor 
intensity and is readily comparable between different samples [33]. The simulations (Figure 2) show that 
NFRET ratios provide measures that are similar to E 32 / E 23 ratios. We have therefore used pixel-by-pixel 

sensitized emission NFRET as the most appropriate parameter to monitor changes in subunit 
stoichiometry [6,18,19,27]. 

Figure 3 presents further analysis of our most complete experiment, previously published as 
supplementary information [27]. For the experiment of Figure 3, when the [32-EGFP mole fraction was 
0.67 and 0.33, the average NFRET was 0.0798 and 0.101, respectively. We tentatively assume that 
these ratios produce nearly pure populations of (a4) 2 ((32) 3 and (a4) 3 ([32) 2 nAChRs, respectively (this 
assumption will be explored further in the next section). These data yield E 32 / E 2 3 =1.27 . Thus, the 

data confirm that the NFRET approach can differentiate between the (a4) 2 ([32)3 and (a4)3([32) 2 receptor 
stoichiometries. Figure 2 A indicates that the data would be satisfied by Ri A ve values between 35 and 
50 A, depending on the assumptions for G and A. Our analyses assume the a4, [32 subunit order 
in which there are two P2<a4 interfaces (Figure 1A). The data would also be consistent with a 
single-interface model (see Supplementary Material). 

Figure 3. Data from NFRET measurements with biased transfection ratios. (A) NFRET 
images for representative cells transfected with 0.5 mole fraction plasmid ratios of 
[32-EGFP and a4-mCherry. The lookup table runs from 0 to 0.2 fractional NFRET. Scale 
bars, 5 um; (B) NFRET data for transfection with 0.67, 0.5, and 0.33 mole fraction plasmid 
ratios of [32-EGFP. For each graph, the heavy line is the measured histogram; this is 
generally obscured by the blue curve which sums the Gaussians fitted to the high and low 
NFRET components (dashed lines). The computed fractional area of high NFRET pixels 
{Whigh) is shown in each case; (C) Average NFRET and Whigh vs. mole fraction of 
[32-EGFP cDNA transfected, from the data in (B). The lefty-axis gives mean NFRET. The 
standard errors of the mean are much smaller than the size of symbols, because each 
sample had > 10 6 pixels. The right j-axis is Whigh', (D) Plots of the CV of the NFRET pixel 
distribution (y-axis) vs. the expression level (x-axis). Each point corresponds to a cell. 
Expression level is quantified as NFRET normalization, the mean of I eGFP sjl mcherrv ) for 

each cell. The fitted curves have the form, y = b(a + yfx) . 
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Figure 3. Cont. 
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2.2.1. Further Analysis Based on Distinct Populations 

We wish to detect changes in the fractions of receptors with each stoichiometry, utilizing the full 
power of analyzing each pixel. The most appropriate way to analyze NFRET distributions uses 
Gaussian components [19,22,30]. We note that two Gaussian components are required to fit the 
pixel-by-pixel distribution of NFRET in most cells. The average NFRET values differ among cells 
transfected in the same dish, and also between subcellular regions of individual cells. These 
characteristics can be observed in Figure 3. At the subcellular level, these variations exceed those 
observed for intracellular images of soluble proteins undergoing FRET [28]. Previous studies suggest 
that these variations arise from the fact that the nAChRs are localized primarily on organelle 
membranes [6]. To characterize the distribution of FRET values among pixels, but without committing 
to a particular source of the variation, we pool the entire distribution of pixels in the entire collection of 
30-55 cells under each condition [6,19,31]. The predicted NFRET 32 / NFRET 2l < 1.5 , so that one does 

not expect to observe two well-resolved peaks; and we found no such case. We therefore fitted pixel- 
by-pixel NFRET histograms to two Gaussian components, with peak (=average) values denoted 

NFRET lugh and NFRET lmv . 

Trends similar to Figure 3A-C were also observed in two less complete experiments. The average 
ratio, NFRET high I NFRET lm . in all experiments was 1.23 and the full range of measured ratios was 1.18 

to 1.25. Thus the NFRET distributions display peaks corresponding to the extremes of NFRET that we 
measured by biasing the cDNA ratios. If the NFRET values for the pure stoichiometries, NFRET^j, and 
NFRET 2 ,i, are greater and less than NFRET/ lig h and NFRET[ 0W , respectively, we would expect the 
NFRET distributions to include peaks or shoulders at these extreme values, outside the range defined 
by the averages for the biased cDNA ratios. However the distributions reveal no such extreme 
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components. These arguments provide confidence that the NFRETh ig h and NFRET[ ow components 
correspond to the (a4) 3 ([32) 2 and (a4) 2 ([32) 3 stoichiometries. 

We calculate the integrals, Inti ow and Inthigh, of the NFRET/ ow and NFREThigh component in each 
case. We then calculate the weighted contribution of the Inthigh component: 



W high = Int high I {Int low + Int high ) 



(10) 



Cells transfected with a mole fraction of 0.5 [32-EGFP cDNA subunits and 0.5 a4-mCherry showed 
a nearly equal proportion of high and low NFRET distributions (Whigh = 0.56). For cells expressing 
0.67 and 0.33 mole fraction of [32-EGFP cDNA, W hig h was 0.16 and >0.98 respectively (Figure 2). 
While there is no straightforward way to provide error estimates for Whigh, note that we require a 
coefficient of determination value of 0.995 to accept the number of components (one or two) in the 
fitted distributions [6,27]. The idea that a pure (a4) 3 ((32) 2 population forms more readily than a pure 
(a4) 2 ((32) 3 population agrees with recent data from stably transfected cell lines [34]. Thus we conclude 
that Whigh increases monotonically, but perhaps not linearly, with the fraction of (a4) 3 ([32) 2 
stoichiometry. Therefore changes in Whigh denote changes in the fraction of a4[32 nAChRs with the 
two stoichiometries. 



2.2.2. Variations in Stoichiometry among Pixels and among Cells 

The formal Gaussian fits in our experimental data (Figures 3 and 4) [6,27] could arise either if 
(a) each component contains a mixture of individual pixels, each containing a pure (a4) 3 (P2) 2 or 
(a4) 2 (P2) 3 population, or (b) individual pixels contain a mixture of (a4) 3 ([32) 2 or (a4) 2 ([32) 3 . The 
subcellular mechanisms that partially segregate these populations are not yet clear. At present our 
ability to distinguish among these is limited by variability in FRET measurements arising in part from 
optical distortions [32]. 

Figure 4. Data from NFRET measurements with pharmacological chaperones. Further 
analysis of a published experiment [27]. (A) NFRET images for representative cells 
transfected with 0.5 mole fraction plasmid ratios of [32-EGFP and a4-mCherry. Incubation 
for 48 h in 0.1 [iM nicotine (left), no added drug (center), and 0.1 \xM cytisine (right). The 
lookup table runs from 0 to 0.2 fractional NFRET. Scale bars, 10 j^m; (B) Average NFRET 
and W/ugh for each incubation, from data like those of Figure IB. The left y- axis gives mean 
NFRET. The standard errors of the mean are much smaller than the size of symbols, 
because each sample had >10 6 pixels. The right j-axis is Whigh- 
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NFRET serves well to compare data among cells with varying expression levels. Nonetheless, a cell 
with a higher expression level is expected to provide a less noisy overall estimate of FRET. Suppose 
that the number of nAChR molecules undergoing FRET in a given pixel is a uniformly distributed 
random variable whose mean varies among cells. (A cell is also a collection of compartments that may 
have their own means; NFRET in intracellular compartments has been reported separately [27]). 
Measurements from a uniformly distributed collection of molecules are expected to have standard 
deviation proportional to the square root of the average number of molecules. In this paper, we extend 
the previous analyses by comparing the distributions across cells by plotting the coefficient of 
variation (CV) of the NFRET signal vs. the expression level. We approximated the expression level 
within a cell by taking either its mean net FRET or its mean normalization value (denominator in 
Equation 8); results were similar for the two metrics, and the second is presented in Figure 3D. In the 
data for 0.5 mole fraction [32-EGFP cDNA, we found that the CV varies as the -0.48 power of 
expression level, which is near the expected inverse square -root relation. Figure 4D displays further 
analyses for mole fractions of 0.33, 0.5, and 0.67 [32-EGFP cDNA. We fitted the distributions to the 
sum of a constant term and an inverse square-root term, y = b(a + yfx) . The small a value in the 

middle plot indicates that the CV shows a nearly ideal inverse square -root dependence on the 
expression level. The larger a value in the first and third panels show that CV has an additional 
constant component. The inverse square-root term represented nearly the entire distribution for 0.5 mole 
fraction, but the distribution had a markedly larger constant term for the two biased transfections. 

Thus, the biased expression levels generate a component of CV that does not decrease with 
expression levels. This extra variation has an unknown source. We note that biased transfections are 
expected to produce an excess of unpaired donor or acceptor fluorophores, depending on the subunit 
ratio. Therefore part of the excess variation may arise from pixel to pixel variations in f D or in e, 
a parameter related to variations in JT^ (Supplementary Figure S2). 




2.3. Pharmacological Chaperoning Stabilizes Alternative a4f32 Stoichiometrics 

We have also studied pharmacological manipulations thought to produce changes in the subunit 
stoichiometry [6,18,27]. Our experiments have typically utilized the cDNA ratio of 0.5 (as in the 
experiments of Figure 3, middle panels). Pharmacological chaperoning is expected to produce less 
complete changes in stoichiometry than the biased cDNA ratios. As expected, the pharmacological 
chaperones have generally produced NFRET averages different from the cells treated with no drug 
(Figure 4), but intermediate between the extremes of Figure 3. Nicotine produced a lower NFRET, 
consistent with many previous studies that nicotine stabilizes the (a4) 2 ((32) 3 stoichiometry. On the 
other hand, cytisine increased NFRET [27]. This finding, novel at the time, has since been 
confirmed by single-molecule measurements of nAChRs on the plasma membrane, using zero-mode 
waveguides [27]. Clearly, cytisine stabilizes the (a4)3((32) 2 stoichiometry. 

From these pharmacological experiments, we have analyzed the pixel-by-pixel data [27]. This 
analysis, similar to that of Figure 3, extracted the values of NFRET low , NFRET high for the two 

Gaussian components, which correspond to (a4) 2 (|32)3 and (a4) 3 ([32) 2 nAChRs, respectively. The ratio 
NFRET hi h I NFRET low is usually similar to the values obtained with biased transfections [27]. The 
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values of Wu s h = Inthi g hl{Inti 0W + Inthigh) for the untreated cells were similar to those of Figure 3: for 
nicotine, less than the untreated cells; and for cytisine, more than the untreated cells. That nicotine and 
cytisine respectively decrease and increase Wu g h, compared with untreated cells, was observed in each 
of five independent transfections like that of Figure 4. The values of Wu s h produced by the 
pharmacological chaperones are less extreme than those for biased transfections, again indicating that 
pharmacological chaperoning by nicotine and cytisine partially shifts the population toward 
either stoichiometry. 

2.4. One vs. Two $2 < a4 Subunit Interfaces 

The data do not allow for a choice between subunit orders that involve one vs. two [32 < a4 subunit 
interfaces; these interfaces are also thought to be high-affinity ligand binding sites. The two-interface 
model is supported by two facts about a4[32 nAChRs: the Hill coefficient for the dose-response 
relation exceeds unity; and concatameric subunits function well when there are two [32 < a4 units 
within a pentamer [35]. Nonetheless, other Cys-loop receptors (5-HT3, GABA A , glycine, and GluCl 
channels) may display different subunit orders. All theory and analysis in the paper assumes that only 
two populations are present (those pictured in Figure 1A). Supplementary Materials provides 
simulations for the single-interface model, but again for only two populations. 

2.5. Requirements for the Procedure 

We outline the optimal pre-requisites for the procedures. First, one requires some assurance that the 
fluorophores do not markedly distort function. For the present fluorescent constructs, several studies 
provide this assurance [6,16-18,29,36]. It is now also possible to achieve functional responses when 
the fluorescent proteins are placed at the extracellular C-terminus [37]. 

One requires a modest density of receptors. We estimate that the nAChRs in our studies are present 
at a density of 3-100 urn" , because we have utilized similar cells and microscopes for single-molecule 
fluorescence [38-40]. Much lower densities would provide insufficient signals. Much higher densities 
could produce "stochastic" FRET. 

One would prefer a system to determine the NFRET values of nearly pure stoichiometries. This 
allows one to define the NFRET values corresponding to E i2 , E 23 directly, as in Figure 2C. We have 

accomplished this via transient transfections with biased cDNA ratios. In the absence of this preferred 
route, one can extract the appropriate NFRET values from Gaussian fits to NFRET distributions. As 
shown in Figures 3 and 4, the Gaussian fits produce a NFRET hjgh I NFRET low value corresponding well 

to directly determined NFRET . 
3. Conclusions 

These experiments analyze FRET-based measurements of changed nAChR stoichiometry in live 
cells with good spatial resolution. We emphasize the detection of altered stoichiometry rather 
than absolute values. Nonetheless, it is rather satisfactory that the experimental values of 
NFRET,. h I NFRET low are simulated by plausible parameters for the structure of a4[32 nAChRs. It is 
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also satisfactory that consistent conclusions arise from two different analyses of the data: NFRET , 
and Whigh- The values for net FRET, E i2 l E 23 , follow the same trends (Supplementary Information). 

Summarizing these conclusions: the biased subunit transfections favor the two stoichiometries in the 
expected directions; pharmacological chaperoning by nicotine favors the (a4)2((32) 3 stoichiometry; and 
cytisine favors the (a4)3([32)2 stoichiometry [2,41]. 

The present simulations and experiments incorporate donors or acceptors within all receptor 
subunits. This tactic provides the largest possible fluorescence signals [42], roughly twice as large as in 
our previous quantitative FRET methods, where donor and acceptor fluorophores were included only 
in a4 subunits, or only in [32 subunits [6,18]. 

3.1. Investigations that Might Use the Procedures 

The experiments produce valid results independently of knowledge about Ri iave . This is fortunate, 
because the intracellular domains of Cys-loop receptors have not been characterized structurally, 
despite the excellent progress on structures for the extracellular binding regions and transmembrane 
domains [43—45]. We emphasize that distance measurements are not the focus of this study; rather, we 
describe useful experimentally determined FRET-related measurements that reveal changes in stoichiometry. 

Which range of Ri, ave values would yield useful data? Because the fluorescent proteins are cylinders 
with diameter of 25 A and length 40 A [46], they could physically approach as close as Ri_ ave ~25 A, 
which is the lower limit of our simulations (Figure 2A) but predicts E J2 / E 23 and NFRET32/NFRET23 
values markedly lower than we measure. At the other extreme, Ri Ave might be --40 A greater than the 
intersubunit distance of the residues at the point of insertion; but Ri Mve > -70 A would probably give 
FRET values below the measurable range. 

The Introduction summarizes the various pharmacological and electrophysiological characteristics 
that depend on the subunit stoichiometry. It is also possible to employ the pixel-based resolution in 
order to distinguish between nAChR stoichiometry in distinct organelles, providing a cell biological 
basis for the partially distinct pixel distributions. For instance, we concluded that the GABA transporter 
GAT1 exists in different multimerization states in perinuclear vs. near-plasma membrane regions [19]. 
The resolution is poorer than the 69 nm square represented by a pixel; therefore one should remember 
that adjacent pixels are highly correlated. Because the stoichiometry is determined in the ER, but 
governs exit from the ER as well as trafficking through the secretory pathway to the membrane, the 
procedures developed here will aid in several investigations. We showed that the FRET procedures can 
distinguish stoichiometry changes during chronic exposure to nicotinic ligands [5,6]. Stoichiometry 
may also change with expression of auxiliary subunits or chaperone proteins [12], with mutations 
causing autosomal dominant nocturnal frontal lobe epilepsy [18], and with mutations of ER retention 
and export motifs [6]. Some of these manipulations also affect the efficiency of nAChR exit from the 
ER [6]. The procedures described here should be applied in additional studies. 

3.2. Other Cys-Loop Receptors 

We expect most Cys-loop receptors to be characterized by a geometry factor nearer to 1 .62 and by 
an asymmetry factor nearer to zero, compared with the a4[32 receptors, for the reason that intracellular 
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M3-M4 loop lengths of other Cys-loop receptors occupy a more uniform distribution. Thus the 
analysis presented here may be generally applied. 
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